Novel quantum phases of dipolar Bose gases in optical lattices 



o 
o 

(N 
>v 



S. Yi 1 , T. Li 2 , and C. P. Sun 1 

1 Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, 100080, P. R. China and 
2 Department of Physics, Renmin University of China, Beijing, 100872, P.R.China 

(Dated: February 5, 2008) 

We investigate the quantum phases of polarized dipolar Bosons loaded into a two-dimensional 
square and three-dimensional cubic optical lattices. We show that the long-range and anisotropic 
nature of the dipole-dipole interaction induces a rich variety of quantum phases, including the 
supersolid and striped supersolid phases in 2D lattices, and the layered supersolid phase in 3D 
lattices. 
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The theoretical prediction [l| and experimental real- 
ization of the superfluid (SF) to Mott insulating (MI) 
transition in atomic system have triggered tremendous 
research activities on ultracold atoms in a optical lat- 
tice. The unprecedented experimental control over both 
atomic and optical lattice parameters makes it an ideal 
platform for studying the physics of strongly correlated 
quantum system. The physics of atomic gases in an 
optical lattice is captured by the on-site Bose-Hubbard 
model [l|, however, it has the limitation of lacking off- 
site density-density correlation, which is essential for the 
emergence of density wave (DW) order and more intrigu- 
ingly the supersolid (SS) phase [E ll> IE IE • Recently, it 
was proposed that the off-site density-density couplings 
can be realized by the bosons in higher bands of optical 
lattices [7j, by interaction mediated through fermions in 
a mixture of bosonic and fermionic atoms 0] , and more 
directly, by the long-range dipolc-dipole interaction be- 
tween dipolar bosons [5|. Among these, dipolar inter- 
action has the advantage of high tunability and experi- 
mental accessibility, especially with the recently achieved 
chromium condensate pi and ultracold polar molecules 
in a 3D optical lattice [101 ] . 

So far, ultracold dipolar Bose gases have attracted sig- 
nificant theoretical interest, for the dipole-dipole interac- 
tion gives rise to new phenomena in Bose-Einstein con- 
densates [ll|, HH, [H, [3 and provides new schemes for 
quantum computing [151 ]. In particular, when loaded into 
optical lattice potential, the dipolar interaction may in- 
duce exotic magnetic orderd if the dipole moments are 
free to rotate in space [IE Qj], EH • 

In the present work, we map out the phase diagrams 
of polarized dipolar bosons in 2D square and 3D cubic 
optical lattices. For a 2D lattice on xy plane, the dipole 
moments may point either to z- or y-axis, the former cor- 
responds a 2D isotropic model as the dipolar interaction 
is isotropic on the lattice plane, while the latter is referred 
to as the 2D anisotropic model for the anisotropic dipo- 
lar interaction on xy plane. For 3D lattices, we always 
assume that the dipole moments point to z direction. As 
we shall show below, the long-range and anisotropic na- 
ture of the dipolar interaction give rise to extremely rich 
quantum phases in the lattice systems. In addition to 
the SS and checkerboard phases found in Ref. [5j, the 



striped supersolid (SSS) phases are also observed in 2D 
anisotropic model. More remarkably, by tuning the sign 
of the dipole-dipole interaction to negative, we discover 
the layered supersolid (LSS) phase in 3D lattices. 

In the presence of dipole-dipole interaction, the Hamil- 
tonian for the extended Bose-Hubbard model reads [j| 
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where b\ is the boson creation operator at site i and 
rii = b\bi is the corresponding particle number operator. 
t is the hopping matrix element between nearest neigh- 
bors, Uo > is the on-site Hubbard repulsion due to 
s-wave scattering, and throughout this paper, the chem- 
ical potential /i is fixed at 0.4£/o- 

The last line of Eq. ([I]) describes the dipolar interac- 
tion in the lattice system with coupling parameters being 
expressed as 
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where / is a factor continuously tunable between — i and 
1 via a fast rotating orienting field [lj|. As shown below, 
nontrivial quantum phases appear for negative / values. 
Cdd = d 2 /(47T£ ) or ^o^ 2 /(47r) for, respectively, the elec- 
tric or magnetic dipoles, with d being the dipole moment 
and £o (no) being the vacuum permittivity (permeabil- 

ity). 
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where w(r — r^) is the localized (on site i) Wannier func- 
tion of the lowest energy band and 9 is the angle formed 
by the dipole moment and the vector (r — r'). Since T>ij 
only depends on — rj = a(Z x ic + l y y + l z i) with a be- 
ing lattice constant and l x ,y.z being integers, it can be 
equivalently denoted as i z )- 

In principal, t, Uo, and {U^ d } are all determined 
by the optical lattice parameters However, to ob- 
tain a complete phase diagram, we shall allow them 
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to change independently except that the set of {U^ d } 
have to be calculated consistently for a given lattice ge- 
ometry. For simplicity, we calculate 2?(; x ; H z.) using a 
spherical Gaussian function of width er, w(r — r^) = 

7r -3/4 (T -3/2 e -(r-r 1 ) 2 /(2T 2 ) ; fo r which p (0Q0) vanish ^ _ 

All results presented in this work are obtained with 
^(ixiyh) corresponding to (<r, a) — (1,3). Similar nu- 
merical results are found with other (a, a) combinations. 
For practical purpose, we choose to truncate l x .y,z a c- 
cording to -l max < l XjViZ < l max with Z max = 2. We 
have also performed simulations using l max up to 6 and 3 
for, respectively, 2D and 3D lattices, which yield similar 
qualitative pictures. To simplify the notation, we define 
a dimcnsionless parameter 

7 = /c dd (47rft 2 a sc /M)- 1 

to measure the relative strength of the dipolar coupling, 
where a sc is the s-wave scattering length and M is the 
mass of the atom. It can be easily estimated that 7 ~ 
0.033/ for 52 Cr atom and 7 « 5.68/ for a typical polar 
molecule [2l|. Finally, the mean- field ground state of 
Hamiltonian (fl]) is obtained using Gutzwiller ansatz [l], Q 
on lattices up to 96 x 96 and 24 x 24 x 24 sites with periodic 
boundary conditions. Similar results are expected for 
even large number of lattice sites. The Gutzwiller ansatz 
has been widely used to study the ground state and the 
dynamics of cold bosons in optical lattices [3, H, HH , as it 
can capture qualitatively the quantum phases, especially 
in 3D lattices [22| . 

2D isotropic model. - As shown in Fig. [H only MI 
and SF phases appear when 7 is negative, which is sim- 
ilar to the on-site Bose-Hubbard model with I7I playing 
effectively the role of the chemical potential. The bound- 
ary between two adjacent MI phases is a horizontal line, 
whose position can be determined analytically as follows: 
the energy per site for an MI phase is 



S/U = P (p - l)/2 - np/Uo + ie M p 2 



(4) 



where p is the mean particle density and 7£dd is the 
dipolar interaction energy per site for MI-1 phase. The 
boundary between MI-p and MI-(p + 1) phases then lo- 
cates at 7 = (p/Uq — p)/[edd(2p + 1)]. Furthermore, a 
MI phase is stable only if 



7 £dd + 1/2 > 0, 



(5) 



hence the system collapses for 7 < —0.41, in very good 
agreement with our numerical result. This stability cri- 
terion also applies roughly to the SF phase. 

For positive 7, the SS and DW phases emerge in ad- 
dition to the usual SF and MI phases. In particular, the 
large lobe denoted as DW-1/2 corresponds to the checker- 
board phase, where sites are alternating empty and singly 
occupied. We note that the checkerboard and SF phases 
are not connected directly, in accord with the Landau 
theory of phase transition, namely, phases with distinct 
symmetry breaking patterns can not be connected with 
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FIG. 1: (Color online) Phase diagram in i-7 parameter plane 
for 2D isotropic model obtained on a 12 x 12 lattice. Various 
MI-p phases are distinguished by the occupation number p. 
DW-1/2 denotes the checkerboard phase. In DW-h (DW-I) 
region, there exist many small DW lobes with the mean den- 
sities higher (lower) than 1/2. The phase transitions between 
various insulating phases are of first order, while all other 
transitions are of second order. This observation also applies 
to Figs. [2] and |4](a). 



each other continuously. Both DW-h and -/ regions con- 
tain many small DW lobes with different mean densities, 
as these DW lobes are very sensitive to t, 7, and the size 
of lattice, it is very difficult to map them out in detail. 

From above discussions, it is clear that the repul- 
sive dipolar interaction (7 > 0) with sufficient strength 
breaks the translational symmetry, while the attractive 
one (7 < 0) always preserves it. The symmetry breaking 
and preserving property of the dipolar interactions can 
be understood intuitively as follows: supposing that site i 
is occupied, it is energetically favorable that its neighbor- 
ing sites are equally (less) populated if the dipole-dipole 
interaction is attractive (repulsive). 

2D anisotropic model. - As dipole moments are po- 
larized along y-axis, the dipolar interaction along y (x) 
direction is attractive (repulsive) for positive 7, while for 
negative 7, the opposite is true. Based on the relation 
between dipolar interaction and translational symmetry, 
we expect to find stripe ordering here. The phase dia- 
gram presented in Fig. [2] indeed confirms our conjecture. 
Considering, for example, the insulating phases in Fig. [21 
the particle density along the repulsive direction becomes 
periodically modulated for sufficiently large I7I, while it 
remains to be a constant along the attractive direction. 
The structure of a striped insulating phase is then com- 
pletely specified by the particle densities over a single 
modulation period as [p±, p2, ■ ■ • , p m ]- 
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FIG. 2: (Color online) Same as Fig. [T] except for 2D 
anisotropic model. The DW phases for 7 < are distin- 
guished by their density distributions. For 7 > 0, both DW- 
A and -B regions contain fine structures. The shaded regions 
represent the SSS phases. The inset shows the zoom-in plot 
of DW-A region, the numbers on the left hand side of the 7 
axis mark the boundaries between two DW phases, and the 
density distributions are specified on other side. DW-B has a 
similar structure except that the occupation number on each 
site is increased by 1. 

There exists a subtle difference between the dipolar in- 
teractions corresponding to negative and positive 7: for 
a uniform density distribution, the overall dipolar in- 
teraction is attractive (repulsive) for the positive (neg- 
ative) 7. Consequently, as the dipolar coupling grows, 
the mean density generally increases for positive 7, while 
decreases for negative 7. The DW- [2,0,0] phase, on the 
other hand, is an exception, as it has higher mean den- 
sity compared to DW-[1,0] phase. This can be easily 
understood from its density distribution, for which the 
attractive dipolar interaction within each lattice row is 
maximized, while the repulsive one between rows is min- 
imized. As we further increase I7I, the system collapses 
for roughly 7 > 0.63 and 7 < —1.27. More interest- 
ingly, the shaded regions, to the right of the striped solid 
phases, represent the SSS phase with a density distribu- 
tion resembling their insulting counterparts. Recently, 
the stripe ordered SS phase was also predicted in the p- 
band Bose-Hubbard model of 2D triangular lattices [23| . 

To gain more insight into the stripe ordering, we re- 
exam the Hamilton (|TJ) in the hard-core limit, where it 
maps onto a spin-i XXZ model 

H = -tJ2(4sJ + h.c) + U dd s ! s J - h * E s * ( 6 ) 

with h z = /i — 7€dd/2. Here the ferromagnetic order- 
ing in xy plane corresponds to the superfluidity in Eq. 
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FIG. 3: (Color online) Mean-field phase diagram for spin- 
^ model Eq. (JSJl with nearest neighbor coupling. In MI-0 
phase, lattice sites are not occupied. DW-i (DW-y) denotes 
the striped solid phase with density modulation along x (y) 
direction. The shaded regions represent the SSS phases. 



([T]), and the modulation in sf corresponds to the or- 
dering of the density. For simplicity, we only take into 
account the dipolar interaction between nearest neigh- 
boring (nn) sites. Figure [3] shows the phase diagram 
obtained through numerical optimization. Compared to 
Fig. [2l the SSS regions here are rather small. Since inclu- 
sion of the next-nearest-neighboring (nnn) coupling does 
not help much, we conclude that the soft-core nature of 
Eq. ([T]) should be responsible for the large stable SSS 
phases in Fig. [2j This also explains the absence of SSS 
phase for / > 0. For isotropic model, nnn coupling is 
the key factor for the presence of SSS phase. It not only 
provides frustration necessary for SS order, but is also 
responsible for the stripe ordering at large nnn coupling 
strength [24[. Here, however, due to the anisotropic na- 
ture of the long-range interaction, the nn coupling alone 
is sufficient for generating the SSS phase. 

3D lattices. - In Fig. [4] (a), we plot the phase dia- 
gram for dipolar bosons in a 3D lattice. A special feature 
we note immediately is that the right boundary of MI-1 
phase is a vertical line. This is because our cut-off scheme 
on lx,y,z makes the dipolar interaction energy vanishes 
for a uniform density distribution on a 3D cubic lattice. 
For positive 7, the dipolar interaction is repulsive on the 
lattice planes perpendicular to z-axis, therefore DW or- 
der emerges on each lattice layer for sufficiently large 7, 
in analogy to the 2D isotropic model. However, as the 
dipolar interaction is attractive along z-axis, the density 
distributions on different lattice layers are identical. The 
detailed structure of DW-C region is rather complicated. 
In general, the mean density first decreases with 7, then 
increases slightly until the system collapses at 7 « 0.63. 
The SS phase has a similar density structure to that of 
DW-C phase. 

For negative 7, one can naturally expect the layered 
phases due to the attractive dipolar interaction on lattice 
layer. In other words, for a layered phase, the particle 
density within each layer is a constant while periodically 
modulated along z-axis. In Fig. 0] (a), the layered in- 
sulating phase are specified by the densities on lattice 
layers over a single modulation period. Unlike the 2D 
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FIG. 4: (Color online) (a) Mean-field phase diagram for dipo- 
lar bosons in 3D cubic lattice obtained on a 6 x 6 x 6 lattice. 
The lower panels show the 7 dependence of superfluid param- 
eters on high (b) and low (c) density layers for t/Uo = 0.01 
(solid lines), 0.02 (dashed lines), and 0.03 (dash-dotted lines). 

case, the only density modulation period we found is 3 
even for a lattice whose size is not a multiple of 3. Same 
as the 2D isotropic case, the onset of instability occurs 
at 7 w -0.41. 



The structure of the LSS phase is similar to that of a 
layered solid phase: each modulation period contains a 
layer with higher density and two layers with equal lower 
densities. Figure 2] (b) and (c) show the 7 dependence of 
the SF order parameter |(6)| on, respectively, the high 
and low density layers. In general, SF parameter in- 
creases with (7) on high density layers, while decreases 
on low density layers. In the strong dipolar interaction 
limit, I (b) I high can be two orders of magnitude larger than 
\(b)\ low If we treat the low density layers as charge reser- 
voir unit, then the resulting structure is quite similar to 
that of a high-T c cuprate superconductor. 

In conclusion, we have studied the quantum phases of 
dipolar bosons in optical lattice. We shown that, due to 
the long-range and anisotropic nature of dipole-dipole in- 
teraction and its tunability, the achievable phases in this 
system are extremely rich, which makes dipolar bosons in 
a lattice an ideal candidate for studying the exotic phases 
in strongly correlated quantum system. The experimen- 
tal detection of SSS and LSS phases in chromium atoms 
relies on tuning the scattering length to below 10 Bohr 
radius, while they should be observable in typical polar 
molecule systems. 
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